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Abstract 

Bleher and Kuijlaars recently showed that the eigenvalue correlations from matrix ensem- 
bles with external source can be expressed by means of a kernel built out of special multiple 
orthogonal polynomials. We derive a Christoffel-Darboux formula for this kernel for general 
multiple orthogonal polynomials. In addition, we show that the formula can be written in 
terms of the solution of the Riemann-Hilbert problem for multiple orthogonal polynomials, 
which will be useful for asymptotic analysis. 

1 Introduction 

Multiple orthogonal polynomials are polynomials that satisfy orthogonal conditions with respect 
to a number of weights, or more general with respect to a number of measures. Such polynomials 
were first introduced by Hermite in his proof of the transcendence of e, and were subsequently 
used in number theory and approximation theory, see pQ, |2j, and the references cited therein. 
The motivation for the present work comes from a connection with random matrix theory. In 
the random matrix model considered in the eigenvalue correlations are expressed in terms of 
a kernel built out of multiple orthogonal polynomials with respect to two weights 

w j (x) = e- v W +a '*, J = 1,2, ai ^a 2 . (1.1) 

A Christoffel-Darboux formula was given in |3J which leads to a description of the kernel in 
terms of the Riemann-Hilbert problem for multiple orthogonal polynomials ^2]. It is the aim of 
this paper to extend the Christoffel-Darboux formula to multiple orthogonal polynomials with 
respect to an arbitrary number of weights. We also allow more general weights than those in 

dnj). 

Let m > 2 be an integer, and let w\, W2, ■ ■ ■ , w m be non-negative functions on R such that all 
moments f^° QQ x k Wj(x)dx exist. Let n = (ni,re2, ■ ■ ■ ,n m ) be a vector of non-negative integers. 
The (monic) multiple orthogonal polynomial of type II is a monic polynomial of degree \n\ 
satisfying 

tat = 0,...,,-l, , = l,...,m. (1.2) 
Supported by FWO Research Projects G.0176.02 and G.0455.04 



Here we define, as usual, \ti\ = n\ + 112 + • • • + n m . 

We assume that the system is perfect, i.e., that for every ti € (NU {0}) m , the polynomial Pa 
exists and is unique, see jHJ. This is for example the case when the weights form an Angelesco 
system or an AT system, see e.g. jllj . 

The multiple orthogonal polynomials of type I are polynomials ' for k = 1, . . . , m, where 

(k) 

has degree < n& — 1, such that the function 



Q h (x)=J2 a{ h\x)Mx) (1-3) 



k=l 

satisfies 

for j = 0, . . . , |n 



/ 



x J Qft{x)dx 



(1.4) 

1 for j = \ti\ — 1. 



The polynomials A\ ' exist, are unique, and they have full degree 

deg^ fc) =n k -l, 

since the system is perfect. 

The usual monic orthogonal polynomials P n on the real line with weight function w{x) satisfy 
a three term recurrence relation and this gives rise to the basic Christoffel-Darboux formula 

g i flW (») = ^-L^^-fa) : ^.MftM (L5) 

where 

hj = j Pj(x)x J w{x)dx. 

In order to generalize the formula (|1.5|) to multiple orthogonal polynomials, we consider a 
sequence of multi-indices no, n\, . . . , n n such that for each j = 0, 1, . . . , n, 

\fij\=j, n j+1 >nj, (1.6) 

where the inequality is taken componentwise. This means that we can go from tij to nj + \ by 
increasing one of the components of tij by 1. We view no, ti\, . . . , ti n as a path from no = (the 
all-zero vector) to an arbitrary multi-index ti = ti n . This path will be fixed and all notions are 
related to this fixed path. Given such a path, we define the polynomials Pj and functions Qj 
(with single index) as 

Pj = Pftji Qj = Qfij +1 - (1-7) 

Our aim is to find a simplified expression for the sum 
n-l 

K n (x,y) = J2 P j(x)QM- (1-8) 

j=o 

To do this, we introduce the following notation. We define for every multi-index ti and every 

k = 1, . . . , 771, 

hf = ! P H (x)x nk w k (x)dx. (1.9) 
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(k) 

The numbers are non-zero, since the system is perfect. We also use the standard basis 
vectors 

ej, = (0, . . . , 0, 1, 0, . . . , 0), where 1 is in the A;th position. (1-10) 
Our result is the following. 

Theorem 1.1 Let n G N and let no,ni, . . . ,n n be multi-indices such that (|1.6j) holds. Let Pj 
and Qj be as in (jl.7|) . Then we have if n = n n , 

n— 1 m 

(x-y)J2 Pj(x)QM = Pn{x)Qri{y) - E -^Pn-e^Qn+eM- (1-11) 

j=0 k=l K-e k 

It is easy to see that reduces to the classical Christoffel-Darboux formula (jl.5[) in case 

m = 1. For m = 2 the formula was proven in 

Remark 1.2 It follows from (|1.11|) that the kernel (|1.8|) only depends on the endpoint n of the 
chosen path from to n and not on the particular path itself, since clearly the right-hand side 
of (jl.lljl only depends on n. 

This fact can be deduced from the fact that for any multi-index k and for i ^ j, we have 

Pj;(x)Q U€l (y) + P Ugi (x)Q u , i+gj (y) = P s (x)Q u+gj (y) + P Ug .(x)Q Ugi+g .(y). (1.12) 

The relation (|1.12|) follows easily from Lemma 13.61 below. 

Remark 1.3 In ljQf , Sorokin and Van Iseghem proved a Christoffel-Darboux formula for vector 
polynomials that have matrix orthogonality properties. As a special case this includes the 
multiple orthogonal polynomials of type I and type II, when one of the vector polynomials has 
only one component. In this special case, their Christoffel-Darboux formula comes down to the 
formula 

n— 1 n+m-ln-1 

(x-y)Y t P j (x)Q j (y)=P n (x)Q n -i(y)- Yl Y, c ^Pj{x)Qk{y) (1-13) 

j=0 k=n j=0 

where the constants Cj & are such that 
fc+i 

xP k {x) = J2c jik Pj(x), 

j=0 

see also ([3.5)1 below. In the setting of it holds that Cj t k = if k > j + m + 1, so that 
the right-hand side of (|1.13|) has 1 + ^m(m + 1) terms. Note that in our formula ()1.11|) the 
right-hand side has only 1 + m terms. 

Another Christoffel-Darboux formula for multiple orthogonal polynomials similar to the one 
in ^0] has been given recently in 0. 

Remark 1.4 As mentioned before, the formula ([l.lljl is useful in the theory of random matri- 
ces. Brezin and Hikami j5] studied a random matrix model with external source given by the 
probability measure 



_L e -Tr(V(M)-AM) dM 



(1.14) 



defined on the space of n x n Hermitian matrices M. Here we have that V(x) = \x 2 , A is a 
fixed Hermitian matrix (the external source), and Z n is a normalizing constant. For this case, 
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we can write M = H + A, where H is a random matrix from the Gaussian unitary ensemble 
and A is deterministic. Zinn-Justin ^3] considered the case of an arbitrary polynomial V. 

The fc-point correlation function R^Xi, . . . , A&) of the (random) eigenvalues of a matrix from 
the ensemble ()1.14p can be expressed as a k x k determinant involving a kernel K n (x,y) 

Rk(Xi, ... ,Xk) = det(K n (Xi, Xj))i<i,j<k, (1-15) 

see ^Sl- Suppose that the external source A has m distinct eigenvalues a\, . . . , a m with respec- 
tive multiplicities n\, . . . , n m . Let n = (n\, . . . , n m ). Then it was shown in |3| that the kernel K n 
has the form (jl.8|) built out of the multiple orthogonal polynomials associated with the weights 

Wj {x) = e -<yi*)-°i*) , j = l,...,m. 

The Christoffel-Darboux formula Ijl.lljl gives a compact expression for the kernel. 

There is another expression for the kernel (|1.8j) in terms of the solution of a Riemann-Hilbert 
problem. This will be especially useful for the asymptotic analysis of the matrix model (jl.l4j) . 
We will discuss this in the next section. The proof of Theorem II .11 is presented in Section 3. 



2 Link with the Riemann-Hilbert problem 



Van Assche, Geronimo, and Kuijlaars found a Riemann-Hilbert problem that characterizes 
the multiple orthogonal polynomials. This is an extension of the Riemann-Hilbert problem for 
orthogonal polynomials due to Fokas, Its, and Kitaev [S]. We seek Y : C \ R — > c( m+1 ) x ( m+1 ) 
such that 

1. Y is analytic on C \ R, 

2. for x G R, we have Y + {x) = Y^(x)S(x), where 



S(x) 



1 W\(x) W2(x) 









1 



w m (x) 





(2.1) 



3. as z — > oo, we have that 



Y{z) = [1 + 



z n 
z~ ni 
z~ n2 

















where / denotes the (m + 1) x (m + 1) identity matrix. 
This Riemann-Hilbert problem has a unique solution given by: 



Y{z) 



Pn(z) Rn(z) 

C\Pfi-e{{ Z ) CiR H _ Sl (z) 



C2Pn- 



(-2 



C 2 Rn-e 2 { z ) 



X z ) c m R R _ Sm {z) 



(2.2) 



(2.3) 
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where Pft(z) is the multiple orthogonal polynomial of type II with respect to the weights 
wi, . . . , w m and Rf t = (Rjt i,Rft 2i ■ ■ ■ ■> Rft,m) is the vector containing the Cauchy transforms 



R H,j(z) 



1 

2tH 



Pfj(x)Wj(x) 



x — z 



dx, 



and 



2vri 



h 



j = 1, . . . ,m. 



(2.4) 



Van Assche, Geronimo, and Kuijlaars ^2] also gave a Riemann-Hilbert problem that char- 
acterizes the multiple orthogonal polynomials of type I. Here we seek X : C\R — > c( wl + 1 ) x ( wl + 1 ) 
such that 

1. X is analytic on C \ R, 

2. for x E R, we have X + {x) = X^(x)U(x), where 



U(x) 



1 

-u>i(x) 1 

-w 2 (x) 1 

-w m (x) 









(2.5) 



3. as z — » oo, we have 



JT(*)=(I + 0(- 



2 





z ni 
z n2 













(2.6) 



This Riemann-Hilbert problem also has a unique solution and it is given by 
iQn(x)^ 2nzA H (z) 

^1 / Qn+ei (x) JT^ kiAft +Sl (z) 

h^nf Qn+e 2 (x) £z k 2 A H+S2 (z) 



X(z) 



2ni I Qn+e m ( x ) 



z-x ^m^4n+e m ( z ) 



(2.7) 



where = (A* , A^ , . . . , A™ ) is the vector of multiple orthogonal polynomials of type I with 
respect to w\,... ,w m , Qh(z) = YJk=i A^\x)w k (x) and 

3 = 1, 



If - - h {j) 



, m. 



(2. 



It is now possible to write the kernel K n (x,y) in terms of the the solutions of the two 
Riemann-Hilbert problems, see also 0. First, we observe that X = Y~ l . If we look at the 
j + 1, 1-entry of the product Y~ 1 (y)Y(x) = X t (y)Y(x), where j = 1, . . . , to, then we find by 
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B and P7|l 



km A 



V) 

n+e, 



m , (fc) 



.GO 



Cl-fn— el (**0 



C-mPft—e^ 



fc=l ^ 



(fe) 
n— ej. 



(2.9) 



where in the last step we used the expressions (|2.4|) and (|2.<S[) for the constants Cj and fcj. 
Multiplying l)2.9j) by Wj(y), dividing by 2ni, and summing over j = 1, ... ,m, we obtain the 
right-hand side of Ijl.llj) . Therefore we see that 

^ m 

(x-y)ET n (x,y) = — £ w j {y)[Y- 1 (y)Y(x)] j+1A 



2vri 



[ Wl (y) ••• w OT (y) ]y^(y)y(x) 



1 







(2.10) 



It is clear that the right-hand side of (|2.1()|) is for x = y, which is not obvious at all for the 
right-hand side of 

In [3] the Riemann-Hilbert problem (|2.1|) - ([2.2|) is analyzed in the limit n — > oo for the special 
case of m = 2, n\ = na, and weights 

Wl {x) = e ' n ^ x2 ~ ax \ w 2 (x) = e - n & 2+ax \ 

The corresponding multiple orthogonal polynomials are known as multiple Hermite polynomials 
[21^3- The asymptotic analysis of (|2.1|) ~ (|2.2|) is done by the Deift/Zhou steepest descent method 
for Riemann-Hilbert problems, see [7j and references cited therein. 



3 Proof of Theorem 11.11 

For the proof we are going to extend the path no, ni, 

fin+k ~ = e fc> k = 1, 2, 



, n n by defining 



, TO. 



We will also extend the definition (|1.7|) by putting Pj = Pf L - and Qj-\ 
1, . . . , n + m. 

3.1 Biorthogonality and recurrence relations 

The multiple orthogonal polynomials satisfy a biorthogonality relation. 
Lemma 3.1 We have 

Pk{x)Qj{x)dx = 6j ik , 



(3.1) 

Qn, for j = n + 



where 8j k is the Kronecker delta. 
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Proof. This is immediate from the definitions (|1.7|) . the orthogonality conditions l|1.4|) of the 
function Qj and 1)1.2(1 of the polynomial P k and the fact that P k is a monic polynomial. □ 

Because xP k {x) is a polynomial of degree fe + 1, we can expand xP k (x) as 
fc+i 

xP k (x) = Y J c j>k Pj{x). (3.2) 

The coefficients can be calculated by Lemma l3.1l bv multiplying both sides of (|3.2|) with Qj(x) 
and integrating over the real line. That gives us 



c j,k = J xPk(x)Qj(x)dx. (3-3) 

The coefficients cj k are if j > k + 2. 

Because of (|3.1() we can write yQj(y) with j < n—l as a linear combination of Qq, . . . , Q n +m-i 
and we have by Lemma 13. II 

n+m— 1 

y<3i(y) = XI c j,kQk{y) for j = o, . . . , n - 1. (3.4) 

fc=0 

Using the expansions (|3,2|) and (|3.4|) for xPk(x) and yQj(y) we can write 

n—l n—l n—l 

( x ~ y)^2 p k{x)Qk{y) = ^xP k {x)Q k {y) - ^2Pk(x)yQ k (y) 

k=0 k=0 k=0 

n— 1 fc+1 n—ln+m—1 

= ^2Y1 c hkPj{x)Qk{y) - ^2 c k,j p k{x)Qj{y). 

k=0 j=0 k=0 j=0 

A lot of terms cancel. Since Cj tk = for j > k + 2, and c njn _i = 1, what remains is 

n—l n+m— In— 1 

{x-y)Y J Pk{x)Qk(y)=Pn{x)QH{y)- £ ^c jtk P J (x)Q k (y). (3.5) 

k=0 k=n j=0 

We also used the fact that P n = Pn and Qn-i = Qn- Note that 1)3. 5 Jl corresponds to the 
Christoffel-Darboux formula of [TJU, as mentioned in the introduction. 
In the rest of the proof we are going to show that 

n+m— In— 1 m TjC") 

X ^ Cj . fc P J ( X )Q fc (y) = ^^ J R^ e . fc ( X )Q ?i+e . fc (y) (3.6) 
fe=n i=0 fe=l ^n-4 

so that ()3.5|) then leads to our desired formula (|l.llj) . 

3.2 The vector space generated by the polynomials P^-g^, • • • , Pn-e m 

For fixed y, the right-hand side of (|3.6|) belongs to the vector space spanned by the polynomials 
of Pfi-g 1 , • • • , Pfi-e m ■ in this part of the proof, we characterize this vector space and show that 
the left-hand side of (|3.6|) also belongs to this vector space V. 

Lemma 3.2 The polynomials P^_g 17 . . . , Pn-e m are a basis of the vector space V of all polyno- 
mials 7r of degree < n — 1 satisfying 

/ ir(x)x t Wj(x)dx = 0, i = 0, . . . , nj — 2, j = l,...,m. (3-7) 
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Proof. By the orthogonality properties (|1.2|) of the polynomials Pn-ei f° r i = l,---,m, it 
is obvious that they belong to V. We are first going to show that the polynomials Pft-ei are 
linearly independent. Suppose that 

aiPa- Sl + a 2 Pti~e 2 H r- a m Pn-e m = (3.8) 

for some coefficients aj. Multiplying (|3.8|) with Wj(x)x nj_1 , and integrating over the real line, 

(i) (?) 

we obtain a, /i- - =0. Since we get a,- = for ? = 1, . . . , m, which shows that the 

polynomials are linearly independent. 

Suppose next that it belongs to V. Put 

bj = — ^ — J Tr(x)x ri: >~ 1 Wj(x)dx 
and define the polynomial 7Ti by 

TTi = biPfl^ + b 2 PH~e 2 H H b m P H -e m - (3.9) 

Then tt\ — ir belongs to V and 

(ttx(x) -ir(x))x n i- l Wj(x)dx = 0, j = l,...,m. (3.10) 

This means that tt\ — tt satisfies the conditions 

(tti(x) — vr(x)) x l Wj(x)dx = 0, i = 0, . . . , rij — 1, j = l,...,m. (3-11) 

Because iri — it is a polynomial of degree < n — 1 and the system is perfect, it follows from 
(|3.11|) that 7Ti — 7r = 0. Therefore 7r = 7Ti, and 7r can be written as a linear combination of the 
polynomials P n _ Sl Pn~e m ■ 

The lemma follows. □ 

Lemma 3.3 For every k = n, . . . ,n + m — 1, we have that the polynomial 

n-1 

n k (x) = } j Cj >k P j (x) (3.12) 
j=o 

belongs to the vector space V . 

Proof. Clearly 71"/% is a polynomial of degree n — 1. Using ()3.2j) we see that 

fc+i 

TT k (x) = xP k (x) - ^2c j>k Pj(x). (3.13) 

j=n 

The representation ()3.13|) of Tr k and the orthogonality conditions ()1.2|) show that 7r k satisfies the 
relations (|3.7|) . so that ir k belongs to V by Lemma I5~2l □ 

Because of Lemma 13.31 the left-hand side of (|3.6j) belongs to V for every y, and so by Lemma 
13.21 we can write 

n+ra-ln-l m 

Y, c ^ p o( x )Qk{y) = J2My) p ft-e^) (3.i4) 

k=n j=0 j=l 

for certain functions <ftj(y)- The next lemma gives an expression for <ftj. We use the notation 
s*o = (all-zero vector) and 

j 

Sj = J2^k, j = l,...,m. 

k=i 
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Lemma 3.4 We have for j = 1, . . . , m, 

= E h iis^Qn+*M- (3-15) 

i=l 

Proof. Rewriting the left-hand side of (|3.14|) using (|3.12|) and (|3.13|) we obtain 

m n+m—1 n+m—1 k+1 

E iv) P n-e d (x) = E xP k( x )Qk(y)- ^ C j,k P j{ x )Qk{y) ■ (3.16) 

j=l k=n k=n j=n 

Now multiply (|3.16|) with x nj ~ 1 Wj (x) and integrate with respect to x. Then the left-hand 
side gives 

The second sum in the right-hand side of (|3,16[) gives no contribution to the integral because of 
orthogonality, and the first sum gives 

n+m—1 / \ m / ,. \ 

£ \J Pk{x)x^w.j(x)dxj Q k (y) =Y1[J p n+i~i(x)x^ Wj (x)dxj Q n +i-i{y). (3.18) 
Because of the choice Q3.1|l and the definition (|1.7|) we have 

Pfi+i— 1 = P ti+Si-i) Qn+i—l = Qn+Si 4 = 1,... , Tfl. 

Then we see that the integral in the right-hand side of (|3.18|) is zero if i > j + 1 and otherwise 
it is equal to hjf^_~ . Then (|3.15|) follows. □ 

3.3 Completion of the proof of Theorem 11.11 

In view of (|3,14|) and (|3.15jl it remains to prove that 

h^Qn+g, = £ hf+s^Qn+sM (3-19) 
i=l 

for j = 1, . . . , m, and then (|3.6f) follows. 

To establish (|3.19[) we need some properties of the numbers and relations between Q- 

(i) 

functions with different multi-indices. We already noted that ^ 0. We express the leading 

(i) 

coefficients of the polynomials Ai in terms of these numbers. 

(i) 1 
Lemma 3.5 The leading coefficient of A~{g. is equal to —gj. 

Proof. Because of the orthogonality conditions Q1.2[) and (|1 .41) we have that 

P n{x)Qn+Sj (x)dx 

P n{x)A^l^ (x)vjj (x)dx 
leading coefficient of A^-^J J Pfi(x)x nj Wj(x)dx 
leading coefficient of -A-?-. I h# 

Yl-\-€.j J Y\, 

and the lemma follows. □ 
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Lemma 3.6 Let j 7^ k. Then we have for every multi-index n that 

Pn{x) = — {Pft+gj — PfL+e k ) = 7^ {Pn+ej ~ Pn+g k ) (3.20) 

n+Sj n+e k 

and 



Qri (fcs {Qn—Sj Qn—Sk) (a\ \Qrl—gj Qn—gu) ■ (3-21) 

h-. - h i - 

n—e k n ~ e j 



Proof. We know that P^ is a polynomial of degree \n\ that satisfies the orthogonality conditions 
(|1.2|) , It is easy to see that Pfi+g- ~ Pri+g k is a polynomial of degree |n| that satisfies these same 
conditions. Because the system is perfect, we then have that 

1 P H (x) = P n+gj (x) - P H+ei (x), (3.22) 

for some 7 £ 1. Multiplying (j3,22|) with x nk Wk(x) and integrating over the real line, we find 
that 

1 n n+ej n+ej 

This proves the first equality of (|3.20|) . The second equality follows by interchanging j and k. 



Next we show (|3.21j) . It is easy to see that Qn-gj ~ Qn-g k satisfies the same orthogonality 
conditions (jl.4|) as Q^. Since the degrees of the polynomials A^_~. — ^_g> fc do not exceed the 
degrees of A~ for i = 1, . . . , m, it follows that 

iQri = Qft—gj ~ Qn—g k i (3.23) 

for some 7 6 R. To compute 7, we are going to compare the leading coefficients of the polyno- 
mials that come with Wk(x). Using Lemma 13.51 we find that 

1 _ 1 1 



h^K - hl k K 

n—e k n—ej—e k n—ej—e k 

This proves the first equality of (|3.2H) . The second equality follows by interchanging j and k. □ 

Now we are ready to complete the proof of Theorem ll.il 

Proof of Theorem 11.11 

In view of what was said before, it suffices to prove (|3.19j) . Fix j = 1, . . . , m. We are going 
to prove by induction that for k = 0, . . . , j — 1, 

k 

hfQn+g, = £ h^Qn+sM + h% k Qn+s k+ g r (3.24) 

i=l 

For k = 0, the sum in the right-hand side of (|3.24l) is an empty sum, and then the equality ()3.24|) 
is clear. 

Suppose that (|3.24|) holds for some k < j — 2. Taking Q3.21JI with n + Sk+i + &j instead of n 
and k + 1 instead of k, we get 

h nls 

Qn+s k+1 +gj = ~~^J) {Qn+Sk+l ~ Qn+s k +gj) ■ 

n+s k+1 
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Thus 

hi+s k Qfi+Sk+ej = ^f!+s k Qn+s k+1 + 4+s i+1 ^+4+i+e, (3.25) 

and using the induction hypothesis (|3,24|) we obtain (|3,24j) with k replaced by k + 1. 

So (B^H) holds for every k = 0, 1, . . . , j - 1. Taking fc = j - 1 in l|33i|) . we obtain (|3~T9"|) and 
this completes the proof of Theorem ll.il □ 
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